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Abstract 

In this paper, we use the generalized Hurst exponent approach to study the multi- 
scaling behavior of different financial time series. We show that this approach is 
robust and powerful in detecting different types of multiscaling. We observe a 
puzzling phenomenon where an apparent increase in multifractality is measured 
in time series generated from shuffled returns, where all time-correlations are 
destroyed, while the return distributions are conserved. This effect is robust 
and it is reproduced in several real financial data including stock market indices, 
exchange rates and interest rates. In order to understand the origin of this effect 
we investigate different simulated time series by means of the Markov switching 
multifractal (MSM) model, autoregressive fractionally integrated moving average 
(ARFIMA) processes with stable innovations, fractional Brownian motion and 
Levy flights. Overall we conclude that the multifractality observed in financial 
time series is mainly a consequence of the characteristic fat-tailed distribution 
of the returns and time-correlations have the effect to decrease the measured 
multifractality. 



1. Introduction 

A traditional assumption, used in the early studies of financial time series, 
considered that returns are independent, Gaussian random variables. However, 
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uncountable number of empirical studies, initiated by Mandelbrot (1963), have 
shown that empirical returns reveal instead very rich and non trivial statistical 
features, such as fat tails, volatility clustering and multiscaling. Prom that times, 
several models have been proposed to mimic the multiscaling behavior of stock 
market returns. For instance, Benoit Mandelbrot, together with his students 
Luarent Calvet and Adlai Fisher, introduced a stochastic process as generating 
mechanism of stock market returns with a multifractal cascade (Mandelbrot et al. , 
1997; Mandelbrot, 1999; Calvet and Fisher, 2002). Such multifractal processes 
provide us with a new model with attractive stochastic properties, which can 
reproduce some stylized facts of financial markets: fat tails, volatility clustering, 
long-term dependence and multi-scaling. However, the practical applicability of 
earlier versions of multifractal models suffers from its combinatorial nature and 
from its non-stationarity due to the restriction to a bounded interval. The most 
attractive feature of these processes is their ability to generate several degrees 
of long memory in different powers of returns. More recently Calvet and Fisher 
( 2002 ) and Calvet and Fisher ( 2004 ) introduced a new family of iterative multi- 
fractal models: the Markov-switching multifractal (MSM) model which preserves 
the hierarchical, multiplicative structure of the earlier models, but possesses ap- 
pealing asymptotic properties. 

In the recent years, there has been an increasing interest in the application 
of the scaling concept to financial markets Miiller et al. (1990); Lux and March- 
esi (1999); Dacorogna et al. (2001); Lux (2004); Carbone et al. (2004). Scaling 
properties in time series have been studied in the literature by using several tech- 
niques. For the interested reader let us briefly mention here some of them such 
as the seminal work of Hurst (1951) on rescaled range statistical analysis and 
the modified rescaled range analysis of Lo (1991), the multi-affine analysis (Peng 
et al. , 1994), the detrended fluctuation analysis (Ausloos, 2000; Bartolozzi et al., 
2007; Di Matteo, 2007), or its generalization, multifractal detrended fluctuation 
analysis (Kantelhardt et al., 2002). The challenge for empirical and theoret- 
ical researches lies in uncovering what scaling laws tell us about the underling 
mechanisms that generate the data. Furthermore, the empirical scaling evidences 
should be used as stylized facts that any theoretical model should also reproduce. 

In addition to this findings, Schmitt et al. (1999) show that the additive 
models like Brownian, fractional Brownian, Levy, Truncated Levy and fractional 
Levy models are not compatible with the properties of financial data and they 
propose the multifractal framework as an alternative. Bianchi and Pianese (2007) 
argue that partition function of generally non-multifractal processes fitted to the 
financial time series behaves as those of a genuine multifractal process. Jiang 
and Zhou (2008) find that scaling behavior of the original financial datasets can 
not be distinguished from those of shuffled time series. Zhou (2009) investigates 
the components of the empirical multifractality of financial returns and finds 
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that temporal structure has minor impact on the multifractal spectrum. More 
recently, Schmitt et al. (2011) find that introduction of Euro had no influence on 
the statistical properties of the fluctuations of the Euro- Yuan exchange rate. 

In this paper, we analyze the multi-scaling properties of different time series 
by means of the generalized Hurst exponent (GHE) which provides a robust 
estimator to compute these scaling properties (Di Matteo et al. , 2003, 2005). 
There are two types of scaling behaviors studied in the finance literature: the 
behavior of the returns distribution tails as a function of the movement size, but 
keeping the time interval of the returns constant; the behavior of some forms 
of volatility measure as a function of the time interval on which the returns 
are measured. In this study we investigate the link between the two in real 
and simulated data series. Furthermore, to distinguish between the effects on 
multifractality from time-correlations and from fat-tailed return distributions we 
apply the GHE on shuffled data series where the time history is destroyed but 
the return distribution is maintained. 

The main part of this paper concerns the study of the source of the multifrac- 
tality in financial datasets. Kantelhardt et al. (2002) point out that in general, 
we can find two types of multifractality in the time series: (i) Multifractality due 
to a broad probability density function; (ii) Multifractality due to different long- 
range correlations of the small and large fluctuations in time. In the first case, 
multifractality can not be removed by shuffling the series. In the second case, 
the corresponding shuffled series should exhibit uni-fractality, since all long-range 
correlations are destroyed by shuffling. In case that both types are present in the 
data, shuffled data should show different multifractality than the original series. 

We contribute to the debate about scaling properties of the financial returns 
with a rigorous statistical analysis of the problem. In particular, we investigate 
the two types of multifractality both on real financial data and MSM simulated 
time series. To test the robustness of our findings, we also compare the results to 
the simulations from a-stable distribution, fractional Brownian motion and au- 
toregressive fractionally integrated moving average model with stable innovations 
which allows us to study the impact of short memory in the heavy tailed process 
with long range dependence. Our study is structured as follows. Sections 2 and 
3 review the Markov-switching multifractal (MSM) and the generalized Hurts 
exponent (GHE) methods. Section 4 reports the empirical and simulation-based 
results describing the two types of multifractality in the data. In Section 5 we 
check the results for robustness by comparing the simulations from a-stable distri- 
bution, fractional Brownian motion and fractional autoregressive moving average 
model with stable innovations. Finally, results are followed by the conclusions 
given in Section 6. 
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2. Markov-switching multifractal model 

In the Markov-switching multifractal model (Calvet and Fisher, 2004; Liu 
et al., 2008; Lux, 2008; Liu et al., 2007, 2008) asset returns are modeled as: 

r t = <?tu t (1) 

with innovations drown from a normal distribution with average zero and unitary 
standard deviation (ut ~ N(Q, 1)) and instantaneous volatility determined by the 
product of k volatility components, or multipliers, M^ 2 \ . . . , and a 

constant scale factor a: 

*? = * 2 ilM?\ (2) 

i=i 

(i) 

where M t is a random variable drawn from a binomial distribution, which is 
characterized by random draws taking two discrete values with equal probability, 
i.e., mo and mi, (with mi = 2 — mo, and 1 < mo < 2). Each volatility component 

(i) 

Ml is renewed at time t with probability 73 depending on its rank i within the 
hierarchy of multipliers and it remains unchanged with probability 1 — 73. The 
transition probabilities are: 

73 = l-(l-7 fe ) , i = l,...,k, (3) 

with 7fc G [0, 1] and b G (1, 00). Different specifications of Eq. 3 can be arbitrarily 
imposed, i.e. as in Lux (2008). By fixing 6 = 2 and 7*. = 0.5, we arrive to a 
parsimonious specification: 

7i = l-0.5 2< ~\ i = l,...,k. (4) 

The dynamics resulting from Eq. 1 is a particular version of a stochastic 
volatility model. With its rather parsimonious approach, this multi-fractal pro- 
cess preserves the hierarchical structure of the original multi-fractal model of 
asset returns while dispensing with its restriction to a bounded interval. 

3. Generalized Hurst exponent 

The generalized Hurst exponent (GHE) method is used to estimate the scal- 
ing exponents and study the multifractality of the empirical data. The biggest 
advantage of GHE is that it combines the sensitivity to any type of dependence 
in the data to a computationally straightforward and simple algorithm 1 , 



1 The algorithm provided by T. Aste can be freely downloaded at http : //www.mathworks . 
com/matlabcentral/ f ileexchange/30076 

4 



The generalized Hurst exponent method aims to extend the traditional scal- 
ing exponent methodology, providing a natural, unbiased, statistically and com- 
putationally efficient estimator able to capture the scaling features of financial 
fluctuations (Di Matteo et al., 2003, 2005; Di Matteo, 2007; Morales et al., 2012). 
It is essentially a tool to study directly the scaling properties of the data via the 
q— th order moments of the distribution of the increments. Different exponents 
q are associated with different characterizations of the multi-scaling behavior of 
the signal X(t). 

We consider the (/-order moment of the distribution of the increments (Barabasi 
and Vicsek (1991)) of a time series X(t): 

Kq[T) ~ (\X(t)\i) ' (5j 

where r varies in the interval between day and T ma x days. The generalized 
Hurst exponent H(q) is then defined from the scaling behavior of K q (r), which 
can be assumed ((Barabasi and Vicsek (1991))) to follow the relation: 

K q {T)~(r) qH M. (6) 

We would like to remind the reader that this framework is based on the as- 
sumption that the process under study has such a scaling property and the scaling 
property stays unchanged across the observation time window. In practice, finan- 
cial time series show evidence of variation of their statistical properties with time, 
and show dependencies on the observation window. The simplest case might be a 
linear drift r/t added to a stochastic variable X(t) = X(t)+r/t with X(t) satisfying 
Eq. 6, To apply scaling analysis to X(t) process, we need to subtract the drift rjt 
from the variable X(t) where r] can be evaluated from the (X{t+r) —X(t)) = rjr. 
This may be viewed as a first-order approximation, crucial for the correctness of 
the results, and such a linear detrending is performed in our computations. 

For some values of q, the exponents are associated with special features. For 
instance H(l) describes the scaling behavior of the absolute values of the in- 
crements. The value of this exponent is expected to be closely related to the 
original Hurst exponent, H, that is indeed associated with the scaling of the ab- 
solute spread in the increments. H(2) is instead associated with the scaling of 
the autocorrelation function. 



4. Distinguishing between two types of multifractality 

4-1. Data description 

We consider daily data for a collection of stock exchange indices: the Dow 
Jones Composite 65 Average Index (Dow) and NIKKEI 225 Average Index 
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(Nik). Foreign exchange rates: British Pound to U.S. Dollar (US/UK), German 
Mark to U.S. Dollar (DM /US) and U.S. treasury bond rates with maturity 1 
year and 2 years, 3 years, 5 years and 10 years (TBI, TB2, TB3, TB5 and TBW 
respectively) all over the period from June 1976 to November 2010. The daily 
prices are denoted as pt, and returns are calculated as rt = ln(pt) — ln(^_i) for 
stock indices and foreign exchange rates and as r t = pt — Pt~i for TBI, TB2, 
TB3, TB5 and TBW. 

4-2. MSM estimation of the parameters 

As first step, the parameters mo and a in the MSM model are estimated, for 
the following hypothetical number of cascades in volatility k = {5,10,15,20}, 
by following the GMM (Generalized Method of Moments) approach proposed by 
Lux (2008) using the same analytical moments as in his paper. Table 1 reports 
the results, where the numbers within the parentheses are the standard errors. 
We observe that the results for k > 10 are almost identical. In fact, analytical 
moment conditions in Lux (2008) show that higher cascade levels make smaller 
and smaller contributions to the moments so that their numerical values would 
stay almost constant. If one monitors the development of estimated parameters 
with increasing k, one finds strong variations initially with a pronounced decrease 
of the estimates which become slower and slower until, eventually a constant value 
is reached somewhere around k = 10 depending on individual time series. Based 
on these estimated parameters, we have proceeded with an analysis of simulated 
data. 

4.3. MSM simulation 

For each stock market index, for each set of estimated parameters (rho, a) 
and for each k, we have simulated 1000 MSM paths and computed the H(q) with 
q = {1, 2, 3} using GHE for empirical series as well as the simulated ones. First 
we verify that Eq. 6 holds for r G [1, T max ]. Then we compute the values of H(q) 
as averages of the best fit exponents corresponding to different T max £ [5,19]. 
The whole procedure is repeated for different stochastic variables X(t) from Eq. 
5, namely X(t) = p t , X(t) = Ylt'=i l r *'l ano - -^(*) = St'=i r t' representing prices 
and volatility estimates. 

We compute the 99% confidence intervals of all exponents using different 
Tmax values and find that they are stable in their values within a range of 10%. 
Figures 1, 2 and 3 show the K q (r) curves for q = [1, . . . , 3] on empirical data, 
while Figures 4, 5 and 6 show examples of scaling of MSM simulated data on 
the Nikkei estimates 2 , We report all results for different stochastic variables pt, 



2 Similar scaling behaviors have been also found for the other simulated time series based on 
estimated coefficients from Table 1 and for different considered stochastic variables. 
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Table 1: GMM estimations of MSM model for different k. This table shows empirical GMM 
estimates of MSM model with different k and numbers in parenthesis are standard errors; All 



empirical returns are demeaned prior to estimation. 
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1.318 


0.012 
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(0.003) 


(0.091) 


(0.004) 
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0.014 
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(0.003) 
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(0.002) 
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1.257 


0.017 


1.205 


0.015 
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0.014 
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(0.029) 


(0.003) 


(0.031) 


(0.003) 


(0.032) 
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0.010 




(0.033) 


(0.004) 


(0.031) 


(0.003) 


(0.032) 


(0.004) 


(0.032) 


(0.002) 


TBI 


1.630 


0.116 


1.596 


0.118 


1.595 


0.119 


1.595 


0.119 




(0.042) 


(0.010) 


(0.041) 


(0.013) 


(0.042) 


(0.012) 


(0.042) 


(0.013) 


TB2 
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(0.013) 
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^2 \ r t' I an( i S r t' an d we can see that scaling holds for all tested series and thus 
it is appropriate to compute the Hurst exponent from the relation in Eq. 6, 

4-4- Multi-scaling of the data 

When studying the scaling exponents for different q, one can distinguish be- 
tween two kind of processes: (i) a process where H(q) = H, constant, independent 
of q; (ii) a process with H(q) not constant. The first case characterizes uni-scaling 
or uni-fractality and its scaling behavior is determined from a unique constant 
value of the Hurst exponent: H. For instance, this might be the case of the self- 
affine process where qH(q) is linear and fully determined by its H, for example 
Brownian motion. The second case, when H(q) depends on q, the process is com- 
monly called multi-scaling, or multi-fractal and different exponents characterize 
the scaling of different q moments of the distribution. Di Matteo et al. (2005) 
confirm that H(l) ^ H(2) ^ 0.5 across wide variety of markets and instruments, 
thus multi-scaling is a characteristic feature of financial stock markets. 

Figures 7, 8 and 9 show the scaling functions r(q) = qH(q) for both empir- 
ical and MSM simulated series for prices, absolute returns and squared returns 
respectively. Their deviations from linearity confirm the multi-scaling behavior 
of all data. This is a clear signature of deviations from pure Brownian motion 
and other additive or uni-scaling models. 
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Figure 1: (Color online) Empirical K q (r) (and shuffled data in red) as a function of r = 1, . . . , 19 
days on a log-log scale for X(t) — p t . Each curve corresponds to different fixed values of 
q = 1, . . . , 3 from upper to bottom line, while upper bold line, middle bold line and bottom bold 
line correspond q = 1, q = 2 and q = 3 respectively, (a) Dow (b) Nik (c) DM/US (d) US/UK 
(e) TBI (f) TB2 (g) TB3 (h) TB5 (i) TB10. Note that axes are fixed for all plots which make 
them comparable. 
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Figure 2: (Color online) Empirical K q (r) (and shuffled data in red) as a function of r = 1, . . . , 19 
days on a log-log scale for X(t) — ~}2 \ r t'\- Each curve corresponds to different fixed values of 
q = 1, . . . , 3 from upper to bottom line, while upper bold line, middle bold line and bottom bold 
line correspond q = 1, q = 2 and q = 3 respectively, (a) Dow (b) Nik (c) DM/US (d) US/UK 
(e) TBI (f) TB2 (g) TB3 (h) TB5 (i) TB10. Note that axes are fixed for all plots which make 
them comparable. 
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Figure 3: (Color online) Empirical K q (r) (and shuffled data in red) as a function of r = 1, . . . , 19 
days on a log-log scale for X(t) = J^ r t'' Each curve corresponds to different fixed values of 
q = 1, . . . , 3 from upper to bottom line, while upper bold line, middle bold line and bottom bold 
line correspond q = 1, q = 2 and q = 3 respectively, (a) Dow (b) Nik (c) DM/US (d) US/UK 
(e) TBI (f) TB2 (g) TB3 (h) TB5 (i) TB10. Note that axes are fixed for all plots which make 
them comparable. 
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Figure 4: (Color online) Example of K q (r) (and shuffled data in red) on MSM simulated data 
(Nikkei estimates) as a function of r = 1, . . . , 19 days on a log-log scale for X(t) = pt- Each 
curve corresponds to different fixed values of q = 1,...,3 from upper to bottom line, while 
upper red line, middle red line and bottom line correspond q = 1, q = 2 and q — 3 respectively. 
Simulated data for (a) k = 5 (b) fe = 10 (c) fc = 15 (d) k = 20. 
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Figure 5: (Color online) Example of K q (r) (and shuffled data in red) on MSM simulated data 
(Nikkei estimates) as a function of r = 1, . . . , 19 days on a log-log scale for X(t) = |fv|. Each 
curve corresponds to different fixed values of q = 1, . . . , 3 from upper to bottom line, while upper 
bold line, middle bold line and bottom bold line correspond q = 1, q = 2 and q = 3 respectively. 
Simulated data for (a) k = 5(b) k = 10 (c) k = 15 (d) k = 20. 
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Figure 6: (Color online) Example of K q (r) (and shuffled data in red) on MSM simulated data 
(Nikkei estimates) as a function of r = 1, . . . , 19 days on a log-log scale X(t) — X] r t'- Each 
curve corresponds to different fixed values of q = 1, . . . , 3 from upper to bottom line, while upper 
bold line, middle bold line and bottom bold line correspond q = 1, q = 2 and q = 3 respectively. 
Simulated data for (a) k = 5(b) k = 10 (c) k = 15 (d) k = 20. 
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Figure 7: (Color online) Scaling functions r(q) = qH(q) as function of q for empirical, MSM 
simulated data for fe = 5, 10, 15,20 for X(t) = p t . (a) Dow (b) Nik(c) DM/US (d) US/UK (e) 
TBI (f) TB2 (g) TB3 (h) TB5 (i) TB10. 
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Figure 8: (Color online) Scaling functions r(q) = qH(q) as function of q for empirical, MSM 
simulated data for k = 5, 10, 15, 20 for = £ |ry |. (a) Dow (b) Nik ((c) DM/US (d) US/UK 
(e) TBI (f) TB2 (g) TB3 (h) TB5 (i) TB10. 
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4-5. Scaling exponents for prices and volatility 

Once verified that the scaling behavior given in Eq. 6 holds, we can now 
discuss the associated average values of generalized Hurst exponent H(q). Our 
approach is to compare the averages of H(q) of empirical time series with sim- 
ulations based on the estimated parameters of the MSM models. The analysis 
is motivated by our previous results in Liu et al. (2008, 2007, 2008), where we 
found some agreement between MSM simulations and empirical data. We focus 
on the generalized Hurst exponent (GHE) for q = {1,2,3} as for q > 3 scaling 
does not hold any longer. We compute H (q) for the empirical data and compare 
them to 1000 simulated series for each set of estimated parameters for different 
values of k. 

Tables 2, 3 and 4 show results for the different stochastic variables, pt, J2 \ r t'\ 
and Yl r t'- These tables report H(q) of empirical series and averages of H(q) 
from 1000 MSM series simulated using the parameter estimates from the Table 
1. If the MSM model describes the scaling behavior of the empirical data well, 
the identity of H(q) averages from empirical and simulated data should not be 
statistically rejectable. These cases are presented in bold. 

Let us start by focusing on the first half of the results titled "original" in 
these tables. The empirical scaling exponents are varying with different type 
of data, they are all different among each other as well as different from 0.5, 
H(l) / H(2) / H(3) / 0.5. We observe that the exponents H(l) and H(3) of 
simulated series with different k are in most cases slightly different from 0.5, while 
H(2) can not be statistically distinguished from 0.5. The fact that H{2) does 
not change with k and is alway equal to 0.5 might be explained by the fact that 
MSM is reminiscent of the scaling property for the moments of absolute value of 
return (price change) rather than price itself. Therefore, we have performed the 
same study with stochastic variables constructed as absolute returns and squared 
returns. Still, as we can see from the Table 2, MSM model does not capture the 
multi-scaling property of the data in most cases. 

Let us now discuss the results for the exponents associated with the absolute 
returns Yl \ r t>\- We can see in Table 3 that all estimated scaling exponents, H(l), 
H(2) and H(3) are in this case greater than 0.5 for both empirical and simulated 
series. In particular, a jump in the values of the exponents is observed between 
k = 5 and k > 5 for all series. We can see that multi-scaling properties of abso- 
lute returns are captured well by the MSM model which can not be statistically 
distinguished from the empirical data for most of the cases. 

Let us finally discuss the results for the exponents associated with the squared 
returns Y r t' ( see Table 4 ) which are very similar to results from the absolute 
values case. In this case, all the values of scaling exponents are much higher than 
0.5 and scaling exponents of simulated data correspond well to the empirical data 
as they can not be statistically distinguished. 
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In summary these results tell us that the MSM model is able to capture 
the multi-scaling properties of the volatility, but it is not able to capture the 
properties of the stock market returns. 
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4-6. The source of multifractality 

Let us now proceed with the same analysis as in the previous section, but 
this time shuffling the data randomly to better understand the source of multi- 
fractality. One can argue that, if the behavior described in the previous section 
is caused only by the long-range correlations, then the randomly shuffled series 
should show unifractality, and AHshuff = Hshuffi^) — Hshuff(3) = 0, where 
Hghuffio) is the generalized Hurst exponent of the shuffled time series. If the 
multifractality is caused only by the broad distribution, than instead we must 
have AH = AH S huff > 0. 

We follow our analysis by shuffling the time series 33 times 3 to achieve ro- 
bustness of the results. In the shuffling procedure the values of the time series are 
put into random order destroying all correlations. We shuffle empirical as well as 
simulated series, we check the scaling behavior and we report the results for the 
shuffled data in Figures 1, 2, 3, 4, 5, 6, Scaling of shuffled series is provided by the 
red lines (color online) and we can confirm that scaling holds for all considered q 
values thus we can proceed in computing the scaling exponents 4 , 

If we look at Tables 2, 3 and 4 where we report the results for the prices, 
absolute returns and squared returns respectively, we can compare the scaling 
behavior of "original" series with the "shuffled" ones. We note that the empirical 
scaling exponents are varying with different type of data even after shuffling, 
although they are generally closer to 0.5. 

Let us start discussing the results for the prices. After shuffling, empirical data 
show even larger degree of multifractality than original ones as AH < AHshuff- 
H(2) generally can not be distinguished from 0.5 thus long-range correlations 
have been destroyed by shuffling and shuffled series from the MSM model seem to 
capture the scaling behavior of the shuffled empirical data closer than the original 
ones. Still, shuffled MSM data show stronger multifractality than original ones 
in terms of differences of the scaling exponents. 

Results of absolute values of returns in Table 3 show similar behaviors. While 
shuffling destroys strong auto-correlations and H(2) values drop from around 0.7- 
0.8 to 0.5, empirical data again show stronger degree of multifractality. Simulated 
series from MSM model show the same behavior. Although the difference is less 
than 0.1 in many cases (statistically insignificant), in lot of cases the difference is 
larger, even around 0.3. Table 4 shows the results from the same analysis done 
on the squared returns, Yl r t<- These results are again very similar but the degree 
of multifractality of original as well as shuffled series is even larger. 



3 Note that we have tried different numbers of shuffled time series from 2 to 100, but results 
do not change. These results are available upon request from authors. 

4 Similar scaling behaviors have been also found for the other simulated time series based 
on all pairs on estimated coefficients from the Table 1 and for different considered stochastic 
variables. 
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Table 5: Differences AH and AHshuff ° n original and randomly shuffled data with different 
stochastic variables, pt, ~}2 \r t i \ and ^2,r t ,. Emp. refers to the empirical exponent values, k — 5, 
k = 10, k = 15 and k = 20 refer to the mean and standard deviation of the exponent values 
based on the 1000 simulated time series. 

Pt 



Emp. k = 5 k = lO k = 15 k = 20 Emp. fc = 5 fc = 10 k = 15 fc = 20 

Dow 0.059 0.031 0.026 0.027 0.027 0.123 0.053 0.089 0.100 0.078 

Nik 0.040 0.047 0.050 0.051 0.053 0.086 0.088 0.170 0.203 0.145 

UK 0.048 0.015 0.011 0.011 0.011 0.067 0.033 0.031 0.030 0.038 

Aud 0.039 0.041 0.034 0.033 0.034 0.083 0.072 0.098 0.096 0.116 

TBI 0.154 0.094 0.098 0.098 0.099 0.246 0.181 0.277 0.281 0.294 

TB2 0.165 0.121 0.125 0.130 0.129 0.230 0.222 0.301 0.371 0.291 

TB3 0.134 0.091 0.095 0.098 0.098 0.199 0.168 0.320 0.359 0.274 

TB5 0.102 0.103 0.109 0.111 0.109 0.162 0.198 0.299 0.319 0.275 

TB10 0.123 0.114 0.123 0.129 0.127 0.151 0.217 0.273 0.406 0.352 

EM 

AH AHshuff 

Dow -0.009 0.063 0.066 0.072 0.073 0.171 0.070 0.122 0.113 0.092 

Nik -0.008 0.086 0.102 0.108 0.108 0.112 0.096 0.159 0.189 0.134 

UK 0.053 0.037 0.035 0.038 0.038 0.084 0.046 0.047 0.044 0.061 

Aud 0.049 0.077 0.078 0.084 0.084 0.095 0.084 0.108 0.110 0.131 

TBI 0.115 0.138 0.150 0.157 0.157 0.258 0.172 0.252 0.266 0.293 

TB2 0.125 0.161 0.177 0.181 0.181 0.241 0.210 0.269 0.335 0.232 

TB3 0.064 0.134 0.148 0.153 0.154 0.210 0.160 0.278 0.309 0.245 

TB5 0.044 0.145 0.157 0.165 0.164 0.172 0.192 0.272 0.291 0.248 

TB10 0.096 0.156 0.171 0.182 0.180 0.178 0.185 0.262 0.373 0.313 

Eg 

AH Shuff 

Dow 0.261 0.237 0.242 0.251 0.253 0.387 0.251 0.343 0.335 0.291 

Nik 0.084 0.274 0.302 0.310 0.311 0.322 0.289 0.359 0.421 0.371 

UK 0.270 0.189 0.183 0.188 0.189 0.285 0.221 0.228 0.205 0.240 

Aud 0.231 0.261 0.263 0.274 0.269 0.300 0.289 0.315 0.331 0.342 

TBI 0.383 0.349 0.371 0.378 0.379 0.477 0.386 0.471 0.482 0.503 

TB2 0.406 0.375 0.403 0.401 0.404 0.459 0.414 0.467 0.529 0.459 

TB3 0.295 0.343 0.368 0.371 0.374 0.436 0.371 0.484 0.523 0.457 

TB5 0.262 0.357 0.379 0.387 0.384 0.374 0.399 0.487 0.512 0.451 

TB10 0.448 0.369 0.393 0.405 0.403 0.398 0.417 0.474 0.566 0.519 
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4-7. Does all multifractality come from the distribution? 

Our results from the comparison between the scaling exponents of the original 
series and shuffled ones suggest that most of the multifractality is caused solely 
by the broad, fat-tailed, distribution of the returns. Indeed if we look at Table 5 
that summarizes the results and brings the differences of scaling exponents, AH 
and AHsh u ff, we can see that in most of the cases, AH < AHghuff- Although 
most of these cases is within standard errors and AH and AHghuf f can not be 
statistically distinguished, in some cases the difference is statistically significant. 
This suggests that there might be a negative bias in the estimation of the degree 
of multifractality brought by unknown time-correlations. This result is valid for 
prices as well as absolute returns and squared returns for both empirical as well 
as simulated data from MSM. 

In the next section, we use several different processes to try to explain this 
type of behavior and possibly find the explanation which will help us to better 
understand the source of multifractality in the financial markets. 

5. Robustness of the results 

In order to demonstrate that our results are robust, we have generated time 
series from fractional motion, random walk with Levy distributed steps (Levy 
flights) and autoregressive fractionally integrated moving average (ARFIMA) pro- 
cess with stable innovations allowing us to put short memory into the processes. 

The AH for the Brownian motion is equal to zero as it is uni-fractal (H(l) = 
H(2) = H(3) = 0.5) and shuffling does not change this result. The AH of 
fractional Brownian motion will also be zero (H(l) = H(2) = H(3) = H) and the 
shuffling will destroy the long-range dependence in this case, and still AHghuff 
should be zero. On the other hand, when we turn to uncorrelated series with fat- 
tailed return distributions, the situation is changing. Kantelhardt et al. (2002) 
discuss the case of Levy distribution and show that expected value of H(q) for 
q > a is \jq and expected value of H(q) for q < a is 1/a. In Barunik and 
Kristoufek (2010), we use GHE method and show on the large finite sample study 
that this behavior holds. Finally, we borrow ARFIMA process of Kokoszka and 
Taqqu (1996) to combine both long-range dependence and heavy tails into single 
process. Moreover, we would like to show also the impact of the simple short 
memory processes on the multifractality estimation. 

5.1. Comparison to simulations from ct-stable distribution 

Stable distributions form a class of probability laws with appealing theoreti- 
cal properties which describe well-known stylized facts such as skewness, excess 
kurtosis and heavy tails. Such distributions are described by four parameters: 
a, (3, 7, 5, where: a is the characteristic exponent and < a < 2 ; f3 is the 
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Table 6: H(l), H(2), H(3) and AH of simulated random variables from a-stable distribution 
with different a. Values refer to the mean and standard deviation of the exponent values based 
on the 1000 simulated time series. Note that for a = 2, we have Gaussian random variables 
with mean zero and variance one. 



ORIGINAL 






ff(l) 


H{2) 


ff(3) 


AH 


a = 


1.2 


0.811 


(0.059) 


0.499 


(0.006) 


0.333 


(0.007) 


0.478 


a = 


1.4 


0.705 


(0.035) 


0.499 


(0.014) 


0.334 


(0.020) 


0.371 


a = 


1.6 


0.626 


(0.040) 


0.500 


(0.007) 


0.340 


(0.013) 


0.286 


a = 


1.8 


0.554 


(0.017) 


0.499 


(0.009) 


0.363 


(0.026) 


0.191 


a = 


2 


0.499 


(0.007) 


0.499 


(0.007) 


0.499 


(0.008) 


0.000 


SHUFFLED 






ff(l) 


H(2) 


ff(3) 


AH 


a = 


1.2 


0.811 


(0.058) 


0.500 


(0.004) 


0.334 


(0.005) 


0.477 


a = 


1.4 


0.704 


(0.036) 


0.499 


(0.006) 


0.334 


(0.008) 


0.369 


a = 


1.6 


0.624 


(0.036) 


0.499 


(0.006) 


0.340 


(0.012) 


0.285 


a = 


1.8 


0.555 


(0.019) 


0.500 


(0.008) 


0.363 


(0.026) 


0.192 


a = 


2 


0.500 


(0.010) 


0.499 


(0.009) 


0.499 


(0.009) 


0.001 



skewness parameter and — 1 < @ < 1; 7 and 5 are scale and location parame- 
ters, respectively. The stable distributions have two tails that are asymptotically 
power laws. 

The a-stable distribution can be described by a characteristic function: 

/ exp(- 7 "|nr[l+ l /3(tan^)sign(K)(| 7 - U | 1 -"-l)]+ I ^) a^l 
9W \ exp(-7|w|[l-M/3fsign(w)ln( 7 |w|)] +itu) a = l, {1 

which is the inverse Fourier transform of the probability density function (Nolan, 
2003), i.e., <f>{u) = E[exp(iuX)]. In order to simulate random stable variables, 
we use a method of Chambers et al. (1976). For all values of the parameters 
a < 2 and — 1 < f3 < 1, we set the parameters (a, /3, 7, 5) to (a, 0, \/2/2, 0), where 
a = {1.2, 1.4, 1.6, 1.8, 2}. We choose parameters with special case of a = 2 being 
Gaussian distribution, so the lower the a, the fatter the tails other parameters 
being equal. This will allow us to show the direct impact of the fat tails on the 
multi-scaling. 

From Table 6, we can see that theoretical expectations hold. For the example 
when a = 1.6 corresponding to the estimates on empirical data (Barunik and 
Vacha, 2010), we have the expected value of H{\) = 1/a = 0.625, and expected 
value of H(3) = 1/3 ~ 0.33 and thus expected value of AH = AHshuff = 0.292 
approximately, which are all in good agreement with what reported in Table 
6 within the standard errors of GHE estimator. In general, Table 6 confirms 
that these relations are well followed and the Hurst exponents (as well as the 
differences) are the same for original and shuffled data. Thus we can see that 
heavy tails of the distribution can directly cause large difference in the H(q) and 
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Table 7: H(l), H(2), H(3) and AH of simulated random variables from fractional Brownian 
motion fGn(d) with different memory parameter d. Values refer to the mean and standard 
deviation of the exponent values based on the 1000 simulated time series. Note that for memory 
parameter d = 0.5, we have Gaussian random variables with mean zero and variance one. 



ORIGINAL 






H(2) 


ff(3) 


AH 


fGn(0.3) 


0.300 


(0.008) 


0.299 


(0.008) 


0.299 


(0.009) 


0.000 


fGn(0.4) 


0.400 


(0.008) 


0.400 


(0.008) 


0.399 


(0.008) 


0.000 


fGn(0.5) 


0.501 


(0.009) 


0.500 


(0.009) 


0.499 


(0.010) 


0.000 


fGn(0.6) 


0.598 


(0.010) 


0.598 


(0.009) 


0.598 


(0.009) 


0.000 


fGn(0.7) 


0.697 


(0.011) 


0.696 


(0.010) 


0.696 


(0.010) 


0.000 


SHUFFLED 




ff(l) 


H(2) 


H(3) 


AH 


fGn(0.3) 


0.499 


(0.009) 


0.499 


(0.009) 


0.499 


(0.009) 


0.000 


fGn(0.4) 


0.498 


(0.008) 


0.498 


(0.008) 


0.497 


(0.008) 


0.000 


fGn(0.5) 


0.500 


(0.008) 


0.499 


(0.008) 


0.499 


(0.008) 


0.000 


fGn(0.6) 


0.500 


(0.008) 


0.499 


(0.008) 


0.499 


(0.009) 


0.000 


fGn(0.7) 


0.499 


(0.009) 


0.499 


(0.009) 


0.498 


(0.009) 


0.000 



spurious multifractality on the finite sample data and, because the data have no 
correlation structure, this multifractality can not be destroyed by shuffling. 

Figure 10 (b) shows also scaling functions of the simulated a-stable random 
variables with different heavy tails. Gaussian case with a = 2 shows straight line 
as expected, as it is unifractal and H(l) = H(2) = H(3) = 0.5. But heavier 
the tails are (the lower the a is), stronger the multi-scaling is. Again, shuffled 
data shows no statistically significant difference as expected. The case of a = 1.8 
replicates well the behavior of the data studied in the previous sections. 

5.2. Comparison with simulations from fractional Brownian motion 

We here perform the same analysis with the data simulated from fractional 
Brownian motion (Mandelbrot and van Ness, 1968; Beran, 1994) with different 
long range parameters H = {0.3, 0.4, 0.5, 0.6, 0.7}. The AH of fractional Brown- 
ian motion with Hurst exponent equal to H is zero as this is a uni-fractal process 
and H(l) = H{2) = H{3) = H. Shuffling will destroy the long-range dependence 
and the exponents must become H{\) = H{2) = H(3) = 0.5, while AHshuff is 
zero. 

Table 7 summarises the results. We can see that GHE method is able to 
estimate all values of Hurst parameter on simulated series very precisely, and the 
result is as expected. As the fractional Brownian motion is uni-fractal, AH is 
equal to zero as H{\) = H{2) = H(3) = H. When we shuffle the series, all 
long-range dependence is destroyed and H(l) = H(2) = H(3) = 0.5. 
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Figure 10: (Color online) (a) Scaling functions r(q) = qH(q) for simulated data from fractional 
Brownian motion for different values of long range dependence. Shuffled series in red. (b) 
Scaling functions r(q) = qH(q) for simulated data from a-stable distributions. 



Figure 10 (a) shows also scaling functions of the simulated fractional Brownian 
motion with different long memory parameters. All cases show straight line 
consistently with the fact that all the simulated data are uni-fractal. Moreover, 
shuffled series are the same as the Gaussian case with H = 0.5. 

5.3. Comparison with simulations from ARFIMA process with stable innovations 

While previous two sections allowed us to compare our results to the long- 
range and heavy tailed processes separately, we would like to complete the analy- 
sis by comparing the results to the heavy tailed process, which exhibits asymptotic 
self-similarity and/or long-range dependence. Moreover, we would like to study 
how simple short memory dependence influences the multifractality measures. 
For this purpose, we use autoregressive fractionally integrated moving average 
(ARFIMA) process with a-stable innovations studied in Kokoszka and Taqqu 
(1996). 

These models, denoted ARFIMA(p, d, q), p,q G N are extensions to linear 
ARIMA(j>, d, q) models, replacing the integer differencing exponent d with an ar- 
bitrary fractional real number < <i < 1 — 1/a and 1 < a < 2 heavy tails param- 
eter from a-stable innovations. A fractional ARIMA process Y = {Y(k),k G Z} 
with a-stable innovations is defined as the stationary solution to the back-shift 
operator equation 

%(B)Y{k) = O q (B)(I - B)- d Z a (k), k£Z, (8) 

where the innovations Z a are an i.i.d. standard a-stable random variables, 
BY{k) := Y(k - 1) and $ p (z) = 1 - faz - faz 2 - ... - ^ P z p , and @ q (z) = 
1 — 0\z — 02Z 2 — ... — 9 q z q with roots outside the unit disk {z G C, L\z\ < 1}. For 
more details about this process, see Kokoszka and Taqqu (1996). 
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To obtain the realizations from ARFIMA process, we use the algorithm 
proposed by Stoev and Taqqu (2004). We simulate the ARFIMA process for 
all combinations of different tails a = {1.2,1.4,1.6,1.8,2} and long range pa- 
rameter d = {—0.2,-0.1,0,0.1,0.2}. The usual H parameter we work with 
through this text translates into H = d + 1/a in this process. First, we sim- 
ulate ARFIMA(0, d, 0), hence we obtain heavy tailed process which exhibits long 
range dependence. Second, we add simple AR(1) dependence to it and simulate 
ARFIMA(1, d, 0) with AR(1) parameter equal to 0.4. Thus we put so-called short 
memory into the heavy tailed long memory process. 

Tables 8, 9 summarize the results for ARFIMA(0, d, 0) and ARFIMA(1, d, 0) 
respectively. Both tables report estimated H(1),H(2) and H(3) for all combina- 
tions as well as H(1),H(2) and H(3) for its shuffled counterparts together with 
the differences measuring the degree of multifractality. We expect that shuffling 
will destroy any dependence structure in the data. Similarly to previous sections, 
we use 1000 realizations to obtain average values of estimates with standard de- 
viations. 

The results suggest that GHE serves as precise method for estimation also 
in the case of heavy tail process exhibiting long-range dependence. When we 
shuffle the ARFIMA(0, d, 0) simulated time series, all long-range dependence is 
destroyed and H(q) estimates are the same as those of pure a-stable process 
reported in previous section. The multifractality of this process is solely caused 
by the fat-tailed distribution. This can be seen by comparison of the results to 
previous section, but also that all shuffled series with the same a parameter has 
the same properties. 

In contrast, Table 9 reporting the results from ARFIMA(1, d, 0) simulations 
reveals very interesting findings. Presence of short memory (AR(1)) process 
translates into the upward bias in the H(q) estimates. Moreover, the difference 
AH is much lower in comparison to the same process without the short memory 
(significant in some cases). When we shuffle ARFIMA(1, d, 0), we arrive to the 
same results as in Table 8 and in previous section which considered only pure 
a-stable process. This means that we destroy all the dependence. 

This final exercise brought us very interesting result, that short memory may 
bring significant contributions to the estimated degree of multifractality. 

6. Conclusion 

In this paper, we present and discuss several new results on the multifractal 
behavior of different financial market datasets. By computing and comparing the 
scaling exponents for empirical and simulated data, we show that the generalized 
Hurst exponent approach is a powerful tool to detect the scaling property of 
financial markets data for different stochastic variables, and it is also a good tool 
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to test the reliability of models such as MSM. By using Monte Carlo simulations 
as well as empirical analysis we show that the MSM model is able to capture 
the multi-scaling behavior of the volatility series, but it is not able to replicate 
the scaling properties of the prices. Large simulations show some agreement of 
the MSM model with empirical data, but most of the multifractality seem to 
be originated from the broad distribution of the series and not from the time- 
dependencies in the model. However, we observed a consistent increment of the 
multifractality after shuffling, i.e. when all temporal correlations are removed 
while the return distribution is preserved. We demonstrated that in some cases 
the increment is significantly larger than the confidence interval. We argue that 
short memory time-correlation in the data may be the cause this increase in the 
multifractality. 

We further compare our results to the simulated series from uni-fractal frac- 
tional Brownian motion containing only long-range dependence, as well as heavy- 
tailed a-stable Levy distribution and autoregressive fractionally integrated mov- 
ing average process with stable innovations and the outcomes confirm our result. 

In conclusion, from the present analysis we can argue that most of the mul- 
tifractality observed in the stock markets data is the consequence to a broad 
fat-tailed distribution of the returns, while some part may be affected by the 
presence of short memory process in the data. Our results are statistically robust 
and hold to various types of financial data; stock market indices, exchange rates 
as well as interest rates. 
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